Where we observe that hehcal turbulence prevails over inertial waves 
in forced rotating flows at high Reynolds and low Rossby numbers 
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We present a study of spectral laws for helical turbulence in the presence of solid body rotation up 
to Reynolds numbers Re ~ 1 x 10^ and down to Rossby numbers Ro ~ 3 x 10"'' . The forcing function 
is a fully helical flow that can also be viewed as mimicking the effect of atmospheric convective 
motions. We test in the helical case variants of a model developed previously (Baerenzung et al. 
2008a) against direct numerical simulations (DNS), using data from a run on a grid of 1536^ points; 
we also contrast its efficiency against a spectral Large Eddy Simulation (LES) (ChoUet and Lesieur 
1981) as well as an under-resolved DNS. The model including the contribution of helicity to the 
spectral eddy dissipation and eddy noise behaves best, allowing to recover statistical features of the 
flow. An exploration of parameter space is then performed beyond what is feasible today using 
DNS. At fixed Reynolds number, lowering the Rossby number leads to a regime of wave-mediated 
inertial helicity cascade to small scales. However, at fixed Rossby number, increasing the Reynolds 
number leads the system to be dominated by turbulent energy exchanges where the role of inertial 
waves is to weaken the direct cascade of energy while strengthening the large scales. We find that a 
useful parameter for partitioning the data is Nc = ReRo — Urms/W^], with Urms, v and Q. the rms 
velocity, the viscosity and the rotation rate respectively. The parameter that determines how much 
the energy cascade is direct or inverse-in which case the cascade to small scales is predominantly 
that of helicity-is linked to Ro. 

PACS numbers: 47.32.Ef, 47.27.Gs, 47.27.Jv 



I. INTRODUCTION 

Turbulent flows at high Reynolds number prevail in 
the atmosphere and oceans, as well as in astrophysical 
settings, e.g. in solar-terrestrial (space weather) interac- 
tions. Because of the limited power of computers, even 
with the petascale efl'orts presently under way, model- 
ing of such flows is needed in order to take into account 
the effect that the unresolved small scales, the so-called 
subgrid scales, have on the large scales, as done for exam- 
ple in the case of numerical weather prediction (Palmer 
2001). Numerous models have been put forward and 
tested over the years, with, as their principal ingredients, 
an eddy viscosity representing the dissipation of energy 
linked to the unresolved small-scale eddies (Smagorinsky 
1963, Kraichnan 1976), and an eddy noise (Leith 1971) 
which mimics the stochasticity of the small scales (see 
Meneveau and Katz 2000 for a recent review). 

In the presence of helicity (velocity-vorticity correla- 
tions) nonlinear terms are weakened. This was invoked 
by Lilly (1986) (see also Anthes 1982) to explain the per- 
sistence of strong storms, as observations of helical fea- 
tures are common in atmospheric flows, e.g. when ana- 
lyzing VORTEX data (Veriflcation of the Origin of Rota- 
tion in Tornadoes Experiment, Markowski et al. 1998). 
Tornadoes can be encountered in two types of struc- 
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tures: one that is laminar (like water spouts) and one 
that is strongly turbulent, typical of storms in the central 
USA. To study the latter, two-dimensional incompress- 
ible studies of axi-symmetric tornado-like vortices in the 
presence of vertical forcing were performed in cylindri- 
cal coordinates with resolutions of 64^ points using eddy 
viscosity (Nolan and Farrell 1999).; an analysis of 120 
simulations in terms of the Reynolds (Re) and Rossby 
(i?o) numbers indicated that the determinant parameter 
was Re/Ro = ^l? jv^, with O the imposed rotation, 
the eddy viscosity and L a characteristic scale. In such a 
study, the rotation is imposed, whereas in a tornado the 
rotation has a somewhat uncertain origin (Rotunno 1984, 
Rotunno and Klcmp 1985): it is thought to be linked to 
the pre-cxistcncc of a downdraft which, together with 
precipitation introduced as in Markowski et al. (2003), 
transports angular momentum to the ground where the 
circulation is then closed (see Wicker and Wilhelmson 
1995, and references therein). More recent studies of tor- 
nado genesis have included the effects of moisture and 
buoyancy, together with stratiflcation and latent-heat re- 
lease; this leads to a more realistic modeling of supercell 
storms and of their maximal vorticity. 

A variety of physical phenomena interact to create the 
dynamics of such flow structures, among which are strong 
local rotation (induced for example by rain in down- 
drafts), strong winds with accelerations several times 
that of gravity, stratiflcation, and boundary layer ef- 
fects. Many experimental, observational, phenomenolog- 
ical, theoretical, and numerical studies have been devoted 
to these interactions and, as the power of instrumenta- 
tion increases, progress will continue to be made. Cru- 
cially, such flows and flow structures are embedded in an 
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atmospheric circulation at very high Reynolds number, 
leading to the formation of a myriad of small-scale intense 
structures and to a turbulent, unpredictable background 
flow. As usual, two approaches to studying these flows 
can be contrasted: on the one hand, models as complete 
as possible must be devised and analyzed in a parametric 
fashion. On the other hand, a reductionist approach, as 
in applied mathematics, calls for drastic simplifications, 
as in the case, for example, of the studies that considered 
the tornado as an axi-symmetric two-dimensional feature 
which were mentioned earlier. 

Here, we propose to reduce the problem in a differ- 
ent way by considering the related issue of hydrody- 
namic rotating turbulence, ignoring the thermodynamics, 
but considering highly resolved three-dimensional fea- 
tures that can help determine turbulent properties that 
may be relevant for complex atmospheric flows, and in- 
sisting on the importance of attaining and modeling as 
high a Reynolds number as possible with, at the same 
time, realistic Rossby numbers. In Mininni and Pou- 
quet (2009bc), using a direct numerical simulation (DNS 
hereafter) on 1536^ points of rotating turbulence forced 
with a helical flow, it was shown that the spectral indices 
for the energy and helicity spectra differ from the clas- 
sical Kolmogorov law. The purpose of this paper is to 
reach higher Reynolds numbers than what was attained 
in this massive DNS, as well as lower Rossby numbers. 
Our motivation is to understand the effect of rotation 
on the flow, the impact of Rossby waves on scaling laws, 
and to quantify the effect of helicity. To that effect, we 
shall first assess the validity of several models against 
this DNS, based on previous work for non-rotating tur- 
bulence (Baerenzung et al. 2008a; see also Baerenzung et 
al. 2008b for the rotating non helical case, and Baeren- 
zung et al. 2009 for the non-helical case of coupling to a 
magnetic field). Armed with this model, we shall exam- 
ine the variation of spectral indices for the energy and 
helicity cascade to small scales as a function of both the 
Reynolds and Rossby numbers, and using a parametric 
study, we shall identify different behaviors depending on 
the dimensionless parameters of the problem. The next 
section gives the basic equations, and tjllll analyzes data 
on temporal evolution and spectra stemming from the 
models and the DNS. Then, mV\ presents the parametric 
study, and finally [jV] pi'csents our conclusions. 

II. EQUATIONS AND MODELS 

The incompressible (V-v = 0) Navier-Stokes equations 
for a flow corresponding to dry dynamics, with velocity 
V and constant density po = 1 read, 

9v 

— +v- Vv + 2f7 X V = -VP + i^Av-l-F^ , (1) 

where P is the total pressure (including the centrifugal 
force), rj is the imposed solid body rotation, taken to be 
in the z direction, and F'^ is a forcing term mimicking 



the input of energy (and helicity) to the system through, 
e.g., buoyancy forces. In the absence of viscosity and 
forcing {v = 0, F"^ = 0), the energy E — (v^/2) and 
the helicity H = {v ■uj/2) (Moffatt 1969) are conserved 
(with a; = V X V the flow vorticity). Conservation laws 
are thought to play an important role in the dynamics 
of turbulent flows. The role of helicity in the absence of 
rotation has been studied in a number of papers (see e.g. 
Brissaud et al. 1973, Kraichnan 1973, Andre and Lesieur 
1977, Hunt and Hussein 1991, Moffatt and Tsinobcr 
1992, Holm and Kerr 2002, Kurien et al. 2004, Chen 
et al. 2005, Jacobitz et al. 2008, Krstulovic et al. 2009): 
the energy and helicity spectra E{k, t) and H{k, t) (aver- 
aged over spherical shells of width Afc =1), follow a Kol- 
mogorov law with E{k,t) ^ k^^/^ and H[k,t) ~ k^^/^. 

Using the rms velocity, Urms , the Reynolds and Rossby 
numbers are defined, respectively, from Eq. ([1]) as 

Re = UrmsLalv , Ro = Urms/2Lon , (2) 

where Lq ~ 2it J E{k)k^^dk / J E{k)dk is the isotropic 
integral length scale. Wc can also evaluate a micro 
Rossby number Ro^ = Wrms/Sil where uJrms is the 
rms vorticity, and which measures the strength of lo- 
cal compared to imposed rotation. In all simulations, 
this number is close to unity or larger; here and in IIIH 
V = 1.6 X 10^^. Note that in the presence of strong ro- 
tation, an inverse cascade of energy to large scales is ob- 
served and the total energy is found to grow with time. In 
this case, Re and Ro given in Table U below are assigned 
with the turbulence statistics evaluated at the onset of 
the inverse cascade, T/. 

The computations we analyze are set in a cubic box 
of length 2tt with minimum and maximum wavenumbers 
kmin ~ 1 and kmax = N/S, N being the number of grid 
points in each direction and using the usual 2/3 dealias- 
ing rule. The code is pseudo-spectral with temporal inte- 
gration using a second-order Runge-Kutta scheme. The 
forcing, of amplitude _Fo is given by a Beltrami flow: 

F^'/Fo = [Bcos{kFy) + Csm{kFz)]x + 
[C cos^kpz) + A sin{kpx)] y + 
[Acos{kpx) + B sin{kFy)] z. (3) 

This so-called ABC flow, an eigenfunction of the curl 
with eigenvalue kp injects both energy and helicity in 
the system. Such flows give a periodic checkerboard of 
overturning rolls with jets at the circulation centers and 
were used in Lilly 1986 to mimic flows in rotating thun- 
derstorms. In all runs we take A = 0.9, B = 1, C = 1.1, 
and kp = 7, so as to cover both the direct cascade of en- 
ergy and helicity to small scales and the inverse cascade 
of energy to large scales. 

In what follows, we consider a DNS run that has been 
performed on a grid of = 1536"^ points, the largest 
ever for rotating turbulence to our knowledge, and sev- 
eral Large Eddy Simulation (LES) runs. The LES tests 
presented in the next section were performed on a grid 
of 96'^ points, while the parametric study of ijlVI was 
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done with LES up to 512^ grid points for the highest 
Reynolds number studied. A statisticaUy steady state of 
turbulence with almost no rotation (il = 0.06) was first 
reached in the DNS, after approximately ten turn-over 
times TiqL = Lo/Urrns- Then, rotation was switched on, 
with f2 between 1.75 and 117 for the parametric study, at 
a time taken as the initial condition (defined as t = 0) for 
all the runs analyzed in this paper (DNS and LES). For 
the LES initialization on grids of Nf^ points {Nl < N), 
the 1536'^ data at i = was truncated at kmax — Nl/2 
(note no de-aliasing is needed in an LES) and the com- 
putations lasted roughly another thirty turn-over times. 

Since the runs use a pseudo-spectral scheme with a 
spectral LES, we write the Navier-Stokes equation in 
terms of the velocity in Fourier space: 

with Efizj the usual antisymmetric tensor. t^(k,t) is a 
bilinear operator for energy transfer written as 

tl{k,t) ^ -iP^p{k)k^ i'/3(p,t)v^(q,t), (4) 

p+q=k 

and the projector Pap{k) ~ Sap — kakp/k"^ enforces in- 
compressibility {kaVa = 0); Einstein summation is used. 
From the equation above, one builds the temporal evolu- 
tion for the energy and helicity spectra E{k) and H{k), 
averaged over spherical shells of radius k = |k|: 

{dt + 2i^k'')E{k,t) = TEik,t)+FE{k,t) , (5) 
{dt + 2iyk^)H{k,t) = TH{k,t)+FH{k,t) , (6) 

with TE,H{k,t) the energy and helicity nonlinear trans- 
fers at wavenumber k and FE.H{k,t) the forcing ex- 
pressed in terms of its symmetric (energy) and non- 
mirror symmetric (helicity) parts. The transfer terms 
involve triple correlations between three Fourier modes 
v(k,t), v(p,t) and v(q, <) with p + q = k, expressing 
the fact that the Fourier transform of the nonlinear term 
is a convolution. Hence, one speaks of triadic interac- 
tions between the modes, and detailed energy and he- 
licity conservation occurs for each such triad. This is 
where the classical closure problem of turbulence arises: 
one now needs an equation for triple moments which will 
involve quadratic correlations, and so on. Many closure 
schemes have been devised over the years, one of the 
most successful in spectral space being the Eddy Damped 
Quasi-Normal Markovian model (EDQNM hereafter, sec 
Orszag 1977). Within the framework of the closure, the 
exact transfer terms Te,h are replaced by approxima- 
tions which we denote TE,H{k,t). The EDQNM closure 
was written for helical turbulence in Andre and Lcsieur 
(1977); it leads to a set of intcgro-differential equations 
for the energy and helicity spectra. The EDQNM ap- 
proximation has been used in atmospheric studies, in ap- 
plications to shear, rotating or stratified flows including 
non-Markovianized versions based on extensions of the 
Direct Interaction Approximation (see O'Kane and Fred- 
eriksen 2008 and references therein). 



The EDQNM balance equations for both the energy 
and the helicity can be expressed in terms of emission 
and absorption. For the energy, the former is related to 
eddy noise (an inhomogeneous term involving energy at 
all pairs of modes p, q such that k = p -t- q as stated 
before); the latter, linear in E{k), gives rise to eddy vis- 
cosity, which corresponds to a drain on the energy at 
mode k, although in some cases it can be negative. 

The EDQNM closure in the helical case leads to two 
extra contributions to the transport coefficients as helic- 
ity contributes to eddy viscosity and eddy noise. Since, 
dimensionally, the eddy viscosity, [i't], in the energy equa- 
tion can be written [i>t] = [i^/k], we can write 

dtE{k) 2{u + vt) k^E{k) - 2i>tk^H{k). (7) 

This uses a short-hand and simplified notation to bring 
out the structure of the model (see Baerenzung et al. 
2008a for details), which omits multiplicative coefficients 
as well as both the resolved scale contributions and the 
eddy-noise contributions. While the EDQNM eddy vis- 
cosity Vt{k,t) ~ J fi{k,p,q)E{q)dpdq depends (through 
the function /i) on an integral of the energy spectrum in 
the small scales and represents the drain of energy from 
the resolved scales due to the unresolved sub-grid scales, 
Vtikjt) J f2{k,p,q)H{q)dpdq gives the contribution of 
small-scale helicity through a different function /2 The 
total eddy noise represents the effects of the small-scale 
helicity as well as the effects of small-scale energy on the 
large scales of the flow. 

In a previous paper (Baerenzung et al. 2009), we tested 
for rotating flows the non-helical version of the model, 
with all pseudo-scalar terms (denoted with a tilde) equal 
to zero. However, it was shown recently that when both 
sizable rotation and helicity are present, the small-scale 
cascade is dominated by helicity (Mininni and Pouquet 
2009ab), whereas the energy mostly undergoes an inverse 
cascade, as is well known. The question then arises as to 
whether a helical model is needed for flows where helicity 
plays a central role in the direct cascade dynamics. To 
that effect, we contrast in the next section results stem- 
ming from three models on 96^ grids: the full helical 
model, hereafter labeled LES-PH, the non-helical model, 
hereafter LES-P, and the ChoUet-Lesieur model, in short, 
CL or LES-CL. For completeness, we also compare the 
1536^ DNS data to an under-resolved DNS on a grid of 
160'^ points (which, using the 2/3 dealiasing rule, has 
kmax = 53, shghtly larger than for the LES runs). 

The Chollet-Lesieur (1981) model introduces an eddy 
viscosity of the form, with kcut = N/2 — 1 a cut-off: 

i^CL{k,t) = Ciy+{k,t)^/E{kcut,t)/kcut ■ (8) 

The quantity v^{k^t) = 14- 3.58(fc/fccMt)* is a dimen- 
sionless cusp function; vchik^t) replaces vt in Eq. ([7]) 
(with Vt = 0). The CL model was derived from the 
EDQNM equations using the fact that C \/E{kcut, t)/kcut 
is the asymptotic expression of nonlocal transfer from 
subgrid to resolved scales and assuming a Kolmogorov 
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spectrum extending to infinity. C in our runs was ad- 
justed with the Kolmogorov constant computed from the 
ABC flow resolved by a DNS run using 512^ grid points, 
VIZ. C = 0.14. 



III. TESTING OF THE HELICAL MODEL 
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FIG. 1: Evolution of the total energy (top), helicity (middle) 
and enstrophy (bottom) for several runs: DNS, 1536^ grid 
points (solid line), under-resolved DNS, 160'' points (circles), 
LES-CL (dotted line), LES-PH (dashes) and LES-P (crosses), 
all on 96^ points. The "DNS Red" (for "reduced") corre- 
sponds to the data of the 1536^ DNS filtered down to 96^ 
grid points. As a reference, the energy of the under-resolved 
simulation at t = 20 is i? = 1.39, and at f = 30 it is 1.69. 



We show in Fig. [T] a comparison of the temporal evo- 
lution of the total energy (top), hehcity (middle) and 
enstrophy (w^) (bottom) for five runs which were done 
specifically for comparison and validation purposes of the 
LES models against the high resolution DNS. They cor- 
respond to three 96"^ LES runs (PH, P and CL) and an 



under-resolved DNS with the same Re and Ro as the 
1536^ DNS; the reduced label stands for 1536^ DNS data 
filtered to 96'^ points. In all that follows, the DNS is 
shown with a solid line, the LES with helical coefficients 
with a thick dotted line and that without the helical con- 
tributions with a thinner dash-dotted line, the ChoUet- 
Lesieur model is shown with thin (grey) bars and finally 
the under-resolved run is displayed with circles. The en- 
ergy displays three distinct temporal phases, all well re- 
produced by the models except for the under-resolved 
DNS), although to various degrees of accuracy. At first, 
the energy grows and then decays (up to t = 10), due 
to an initial adjustment as the Coriolis force gets sud- 
denly larger at i = 0. Even though the Coriolis force 
itself does not input energy in the system, inertial wave 
resonances render the flow quasi bi-dimensional and, af- 
ter t — 10, the subsequent growth of the energy is then 
attributed to the onset and further development of an 
inverse cascade of energy with a transfer to scales larger 
than the forcing scale Lp. The transition times between 
these phases are well reproduced by all models, and so 
is the growth rate in the inverse cascade, except for the 
LES-P model. The under-resolved run stands on its own: 
the lack of small-scale dissipation produces an unphysical 
growth of energy and its results, clearly unreliable, will 
not be commented upon further other than to say that 
with insufficient numerical resolution, a model indeed is 
needed to mimic the effects of the unresolved scales. 

The helical model LES-PH is the one closest to the 
DNS, whereas the CL model is too dissipative, because 
it does not take into account the weaker nonlinearities 
due in part to the partial Beltramization of the flow. 
The temporal evolution of the total helicity follows the 
same pattern, except that all the injected helicity under- 
goes a direct cascade as can be verified by computing its 
flux (Mininni and Pouquet 2009b), and therefore its total 
value as a function of time attains a statistically steady 
equilibrium; when examining the errors (not shown) the 
LES-PH model is best overall, and the CL model is worst. 
We also note that, at early time and up to the onset of 
the inverse cascade of energy for t ~ 10, the non-helical 
LES-P model behaves consistently better. This may in- 
dicate that as long as the inverse energy cascade has not 
begun and we are in a regime dominated by waves, the 
direct cascade to small scale is as much energy as it is 
helicity and thus the helicity in the model is not as essen- 
tial. Overall, the normalized error for the total helicity 
using the LES-PH model is between 10"'^ and 10~^, rea- 
sonable values considering the very large difference in the 
linear grid resolution (by a factor of 16). 

The comparisons between models and DNS for the to- 
tal dissipation T>{t) = 2v (uP'{tyj (see Fig. [l]bottom, not- 
ing all runs have the same u) show a large discrepancy 
which can be partially removed when filtering the DNS 
data down to the resolution of the LES runs. Dissipation 
has its largest contribution from small scales which are 
not resolved in the LES runs (see also Fig. [3] below for 
a comparison of dissipation between LES and unfiltered 
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DNS data for different values of v and than those used 
here). When comparing to the filtered DNS displayed 
with triangles in Fig. [U again LES-P and LES-PH give 
the best results, while the CL model over-estimates the 
dissipation by a consistent amount (« 30%). We note 
that, in the case of the Smagorinsky model, a study of 
helical flows (Li ct al. 2006) shows an under-estimation 
of helicity dissipation by 40%. 




FIG. 2: Energy (top) and helicity (bottom) spectra averaged 
during the inverse cascade, from t = 20 to f = 30; same 
notation as in Fig. [1] Insets: zoom on the large scales; note 
the inadequacy of all models in this range of wavenumbers 
except for LES-PH, i.e. when taking into account the helical 
contributions to sub-grid modeling. 

When comparing spectra for energy and helicity, the 
same conclusions arise, with a better behavior for the 
helical model LES-PH. We show in Fig. [2] the spectra 
for energy (top) and helicity (bottom), averaged over an 
interval of time after the start of the inverse cascade of 
energy. All spectra stop at the maximum wavenumber of 
the LES, fc^f^ = 48 and do not show a dissipation range 
by construction (note, however, that in our approach, 
we keep the (bare) viscosity in the equations, that of the 
DNS run on the grid of 1536"^ points). In the small scales, 
the LES-PH model reproduces very closely the DNS up 
to fc^ff , particularly so for the helicity. In the inverse 
cascade of energy, a good agreement is obtained as well, 
with the CL model being less efficient, a feature already 
noted on the temporal evolution of the total energy. The 
helicity does not undergo an inverse cascade for any of 
the runs performed in this paper, but the agreement with 
the DNS spectra is less striking than for energy. 

For scales larger than the forcing scale Lp^ E{k) ap- 



pears to follow a classical Kolmogorov law. There is 
evidence in other work for much steeper power laws in 
the large scales of rotating turbulence, viz. E{k) ^ k^^ 
(Smith and Lee 2005, see also Tran and Bowman 2003), 
when computing in periodic boxes with variable aspect 
ratio and using hyperviscosity. The differences may be 
linked to the effect that wave interactions have on large 
scales on the one hand, and to a subsequent direct cas- 
cade of enstrophy because of the bi-dimensionalization of 
the flow on the other hand (Smith et al. 1999). These 
results, which may also be Rossby dependent (Chen et 
al. 2005) are, however, beyond the scope of the present 
study which focuses on small-scale properties. The spec- 
trum H{k) at large scales is flat and no clear power law 
can be identified. Note that the realizability condition 
\H(k,)\ < kE{k) does not present a dynamical constraint 
for the helicity as fc ^ in the inverse cascade. 

Finally, we note that LES-PH reproduces well the type 
of structures observed in the flow, and in particular the 
spatial juxtaposition of laminar updrafts that are fully 
helical (Beltrami core vortices) and a tangle of small-scale 
vortex filaments ordered in columnar structures (Mininni 
and Pouquet 2009c). The size and shape of eddies in 
the horizontal plane are correctly captured, even though 
phase information is partially lost because of the intrinsic 
stochasticity and lack of predictability of the flow. 



IV. TURBULENCE VERSUS WAVES 

A. The procedure 

Using the LES-PH model tested in Sjlll we now pro- 
ceed to a parametric study of forced helical rotating tur- 
bulence in terms of both the Reynolds and Rossby num- 
bers, examining the scaling laws of energy and helicity 
in the direct cascade. The Reynolds number varies from 
4300 (for smaller values of the Reynolds number it is dif- 
ficult to ascertain the inertia! index of the spectra) to 
1.1 X 10^ , and the Rossby number varies from w 12 to 
0.003. Most runs were done on a grid of 192^ points, al- 
though both lower and higher resolution runs were per- 
formed in order to evaluate to what extent the results 
depended on the resolution of the LES. All runs were 
started in the same fashion, using the same developed 
turbulence state with weak rotation (fi = 0.06) as an 
initial condition, and with the forcing via the ABC flow 
as given in SjlT]with hp ~ 7. Data pertinent to the runs 
can be found in Table [D In the presence of moderate 
to strong rotation, an inverse cascade develops in which 
case the parameters listed in the table were computed 
at the time of the onset of the inverse cascade T/; note 
that only two runs (R3 and R5z) do not show an inverse 
cascade. Note also that LOrms depends on the resolution 
since the enstrophy spectrum peaks in the vicinity of the 
dissipative scale-which is in general not resolved in the 
LES-for energy spectra shallower than fc~^ , leading to an 
under-estimation of the micro Rossby number cjrms/'^^ 
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TABLE I: Parameters of the runs using the LES-PH model; u is the viscosity, SI the imposed rotation, uorms tire rms vorticity 
at Tj, where Tj is the time of onset of the inverse energy cascade; Ro and Re are the Rossby and Reynolds numbers, Urma the 
rms velocity, and Lqjf the integral scales at Ti and t — 30. Data is averaged G [Ti,Ti + 5]. Compared to their values at Tj, 
Urma has increased by roughly 20% at t=29 and u^ma has decreased by between 10% and 20%, except for run R3 for which 
both are stationary, and for run R5F, for which ujrma has decreased by 30%. In the last three columns, the spectral index for 
the energy e, the sum of energy and helicity indices e + /i, and the normalized flux ratio nif/[fcFnE] are given at TV. R3 and 
R5z have no discernible inverse cascade and should have a Kolmogorov spectrum except for the bottleneck that renders spectra 
shallower. The fiducial run R5F is the 1536''^ DNS. Runs of Fi gures 1-2 are not included, except for the DNS run R5F and the 
LES-PH run R5a. 



Run 


Res. 




Q. 




Ti 


10^ i?o 


IQ-^Re 


Urms 


Lo.i 


Lo,F 


e 


e + h 


n„/[kFnE] 


Rl 


192 


25.0 


4.5 


10.5 


17 


4.6 


8.3 


0.9 


2.24 


3.04 


2.1 


4.1 


1.5 


R2 


192 


25.0 


9.0 


11.7 


8 


2.9 


9.3 


1.1 


2.10 


4.30 


2.2 


4.2 


1.7 


R3 


192 


16.0 


1.8 


16.0 


14 


24.9 


4.3 


0.8 


0.88 


0.99 


1.7 


2.8 


1.1 


R4 


192 


16.0 


4.5 


12.4 


14 


6.1 


9.7 


0.9 


1.67 


3.79 


1.9 


3.6 


1.3 


R5F 


1536 


16.0 


9.0 


15.3 


9 


2.9 


14.7 


1.1 


2.14 


4.32 


2.0 


3.9 


1.9 


R5 


192 


16.0 


9.0 


12.8 


8 


3.3 


13.3 


1.1 


1.88 


4.55 


2.1 


3.8 


1.8 


R5a 


96 


16.0 


9.0 


8.5 


10 


2.7 


15.3 


1.1 


2.23 


5.21 


2.2 


3.8 


1.9 


R5z 


192 


16.0 


0.1 


21.8 


7 


1166 


3.6 


0.9 


0.64 


0.64 


1.6 


2.5 


1.1 


R6 


192 


16.0 


18.0 


10.4 


10 


1.6 


16.7 


1.2 


2.16 


4.58 


2.1 


3.9 


2.5 


R7 


192 


16.0 


36.0 


10.8 


10 


0.8 


19.4 


1.3 


2.37 


4.50 


2.1 


3.8 


2.7 


R8 


192 


16.0 


117.0 


12.5 


8 


0.3 


14.9 


1.3 


1.77 


2.99 


2.2 


4.1 


2.3 


R9 


192 


11.9 


18.0 


11.1 


10 


1.4 


25.0 


1.2 


2.41 


4.76 


2.1 


3.9 


2.3 


RIO 


192 


10.2 


42.4 


11.7 


10 


0.6 


32.0 


1.3 


2.46 


3.93 


2.1 


3.8 


2.6 


Rll 


192 


8.0 


9.0 


14.4 


8 


3.8 


24.9 


1.2 


1.72 


3.79 


1.9 


3.3 


1.7 


Rlla 


96 


8.0 


9.0 


7.6 


14 


1.9 


44.9 


1.1 


3.22 


5.87 


2.3 


4.0 


1.9 


R12 


192 


8.0 


18.0 


11.1 


10 


1.3 


43.3 


1.3 


2.75 


4.70 


2.1 


3.8 


2.5 


R13 


192 


8.0 


36.0 


10.9 


10 


0.6 


47.5 


1.3 


2.87 


4.08 


2.0 


3.7 


2.7 


R14 


384 


8.0 


36.0 


15.1 


10 


0.9 


36.6 


1.4 


2.14 


4.18 


2.1 


3.9 


3.2 


R15 


384 


5.3 


72.0 


20.9 


6 


0.6 


42.3 


1.4 


1.58 


3.54 


2.0 


3.6 


3.0 


R16 


192 


5.0 


18.0 


12.4 


9 


1.5 


62.6 


1.3 


2.44 


4.08 


2.0 


3.5 


2.5 


R17 


192 


4.5 


36.0 


15.0 


6 


1.0 


53.8 


1.4 


1.79 


3.91 


2.0 


3.4 


2.8 


R18 


192 


2.5 


36.0 


13.9 


8 


0.9 


108.8 


1.3 


2.02 


3.79 


1.9 


3.2 


3.4 


R19 


512 


1.6 


9.0 


22.7 


8 


3.7 


125.5 


1.2 


1.74 


4.16 


1.8 


3.4 


2.1 



by roughly 50% in some cases. 

We first contrast in Fig. [3] the temporal evolution of 
energy, dissipation and helicity for several runs. Times 
are indicated in units of eddy turnover time, and as the 
rotation rate increases, the inertial wave time decreases 
in the same manner so that the same physical time on 
the plots corresponds to a larger number of wave periods. 
At fixed viscosity (left column), the growth rate of the 
energy associated with the inverse cascade is rather in- 
dependent of rotation except for the lowest rotation rate. 
Otherwise, the time scales are the same for the onset of 
the cascade itself. This is valid as well at fixed rotation 
rate and variable Re as seen in the right column, with 
now comparable energy levels. The discrepancy in the 
dissipation between the DNS and LES at early times has 
already been commented upon in the preceding section 
(see also Fig. [T]). The dynamics of helicity is very similar 
for all these runs, with again variations in the amplitude 
with Re but not so with Ro. 

The main objective of the parametric study is to mea- 
sure the spectral indices of energy and helicity, e and h, 
as well as the helicity to energy flux ratio normalized by 
kp (see below). To measure e, h, a fit is performed in the 



inertial range of E{k) and H{k) of the form: 

E{k) - H{k) - k-'' . (9) 

The spectra for the fit were computed in three ways: (i) 
using the spectra at the time of the onset of the inverse 
cascade T/; (ii) as a time average over an interval At = 5 
starting from t = Ti; and (iii) as an average in the in- 
terval t £ [Ti,2Tj]. Although the actual values of the 
spectral indices were observed to depend slightly on how 
these estimations were performed, it was observed that a 
classification based on whether the sum satisfied e-|-/i ~ 4 
(which will be associated with a wave-dominated regime, 
and corresponds to the black dots in Fig. ^ or e + h ^ 4 
(which will be associated with helical turbulence, and 
corresponds to the gray dots) remains unchanged irre- 
spective of the estimation criteria. The choice finally used 
to define the spectral indices was thus to perform a tem- 
poral average over a not too long time interval (method 
(ii)) to avoid problems associated with non-stationarity 
because of the inverse energy cascade and the energy ac- 
cumulation at large scale. Once the time averaged spec- 
tra were obtained, the least squares fit to get the indices 
was done in the range k £ [15,45]. We also performed 
a spectral fit in the interval k G [19,45] and k G [15,60] 
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FIG. 3: Temporal evolution of energy (top), dissipation (mid- 
dle) and helicity (bottom); the solid line in each plot repre- 
sents the DNS run R5F = 1.6 x lO"*, Q = 9), while all 
other runs are performed with LES-PH. Left column: Fixed 
viscosity {u = 1.6 x 10~*) with Q varying between 4.5 and 18 
by factors of 2 for runs R4 (dash line), R5F, and R6 (dotted 
line). Right column: Fixed rotation rate {Q, = 9) and f de- 
creasing from 2.5 x 10"'' to 8 x 10"^ for runs R2 (dash line), 
R5F, and Rll (dotted line). 



but saw no measurable difference in the results. 

The choice to measure e + h has to do with the two 
expected regimes for helical rotating turbulence. If rota- 
tion is weak, a regime close to Kolmogorov occurs, where 
both the energy and the helicity cascade directly to small 
scales with a k~^/^ spectrum. In this case, which we 
refer to as "helical turbulence", e + h = 10/3. When 
rotation is strong, waves slow-down the energy transfer 
giving E{k) ~ in the non-helical case (DubruUc and 
Valdetarro 1992) and e + /i = 4 in general if helicity is 
present (see Mininni and Pouquet 2009ab). This case, 
we refer to as the "wave-dominated regime" . 

An important ingredient for the helical regime e+h ~ 4 
to occur is that the energy flux to small scales be negligi- 
ble compared to the helicity flux, so that the direct cas- 
cade is dominated by helicity. This can be independently 
verified by computing the energy and helicity fluxes: 
Ue{K) = Ek=,TE{k), and UniK) - ELiTnik). Af- 
ter averaging over the same time interval used to measure 
the indices e and /i, a dimensionless ratio Hh /{kp^E) 



can be computed over the same range of wavenumbers 
k g [15,45]. When examining this ratio function 
of both the Rossby and Reynolds numbers, it is found 
that at low rotation rate the ratio is close to unity, as ex- 
pected, but it is clear that it increases with rotation rate, 
confirming a previous analysis using DNS (Mininni and 
Pouquet 2009a); this ratio increases as well as the tur- 
bulence strengthens (see Table H]). Thus, one can expect 
that helical rotating flows in the wave-dominated regime 
will be characterized by helicity dynamics at small scale. 

Before proceeding to the details of the parametric 
study, we give a few examples of the different behav- 
iors that arise as the Reynolds and Rossby numbers are 
changed. While for most of the flows, the temporal evo- 
lution of global statistics is rather similar (see Fig. [3]), 
the resulting spectra in the runs do differ. This is illus- 
trated by Fig. m which gives (top) the product of spec- 
tra E{k)H{k) for two different runs compensated by the 
two laws discussed above: fc* (solid line with circles) and 
^10/3 ((jash line). Individual energy and helicity spectra 
are also shown in the figure for the same two runs (bot- 
tom), compensated by laws that are the best fit to the 
data. Even though the resolution of the LES is modest, 
the fact that the spectra do not need to display a dis- 
sipation range (dissipation being taken care of through 
the effect of the eddy viscosities) allows for a good de- 
termination of spectral indices. All the spectra shown in 
the figure correspond to a time average from t ~ Tj to 
t = 30. The run on the right of Fig. 2] shows a scal- 
ing which seems to be close to the classical Kolmogorov 
law {e + h ^ 10/3), whereas in the run on the left the 
e + h = A law emerges rather convincingly, with an en- 
ergy spectrum measurably different from a e = 5/3 law, 
the best fit giving E{k) ~ k'^-^, H{k) ~ k'^-^. 

As previously mentioned, spectral indices vary with 
time as seen in Fig. [5l where a fit to obtain e + h was 
performed in the same range of wavenumbers for the in- 
stantaneous spectra as a function of time, using the same 
five runs as in Figure [31 the solid and dash lines repre- 
sent respectively the law derived in Mininni and Pouquet 
(2009a) and the dual Kolmogorov law and are provided 
for reference. Apart from a slower evolution toward an 
established spectral law in the case of lower rotation rate 
(middle row, left-most), the data is compatible with a 
breaking of universality in helical rotating turbulence; 
however, other spectral indices cannot be ruled out (ob- 
serve the middle column). More numerous as well as 
better resolved runs may be needed to assess this point 
further; there may be some ambiguity in assessing inertial 
indices when the runs are close to the transition observed 
in Fig. [6] (see ijlVBI below), and furthermore there could 
be other, as yet unknown, dynamical regimes. 



B. The emergence of two regimes 

With the previously discussed caveats in mind, we now 
examine all runs reported in Table U and classify them 
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FIG. 4: Left column: Compensated spectra for E(k)H{k) 
for run Rl; the bottom plots are the individual spectra E{k) 
(solid) and H(k) (dash), compensated respectively by k^'^ and 
k^'^ . Right column: Same for run R17 (higher rotation and 
Reynolds number), with below the spectra now both compen- 
sated by a Kolmogorov law. 





R2 






R5F 


R6 
















R11 





10 20 30 



FIG. 5: Sum of the energy and helicity indices n{t) = — (e + /i) 
as a function of time. R5F is the large DNS and all other cases 
are LES. Runs 2, 5F and 11 have the same imposed rotation 
(f2 = 9) and decreasing viscosity whereas runs 4, 5F and 6 
have the same viscosity [v — 1.6 x 10"*) and increasing rota- 
tion. The resulting dynamics leads to Reynolds and Rossby 
numbers not quite so well ordered (see Table): for the runs 
in the vertical, the Rossby numbers vary by 50%, and the 
Reynolds numbers by a factor ^ 2.3, whereas for the runs in 
the horizontal, they vary respectively by a factor 6 and 1.8. 
Axes have the same scales in all plots; the values of —4 and 
— 10/3 are shown as solid and dashed horizontal lines. 



based on the resulting sum of spectral indices, e + h. A 
first look at the tabic indicates that, at fixed Reynolds 



number, a transition toward the wave-dominated regime 
takes place as the rotation rate increases, whereas at fixed 
rotation once the eddy turn-over time becomes smaller 
than the inertial wave period, Kolmogorov-like turbu- 
lence takes over. It is not clear from these simulations 
whether the solution e = /i = 5/3is favored, or whether a 
solution like the one postulated in Brissaud et al. (1973) 
for helical turbulence emerges, the differences between 
inertial indices being too minute and the resolution of 
the simulations not large enough. 

In order to check consistency of the results, R13 and 
R14 were run with same [v, 17] but different resolution; 
the corresponding Reynolds and Rossby numbers differ 
respectively by 6% and 15%, whereas the flux ratio differs 
by 18% (see Table U). However, these two runs have the 
same inertial index dynamics. Other convergence tests 
have been performed that are not reported here, leading 
to similar conclusions. 

A demarcation seems therefore clear: at fixed Rossby 
number, a higher Reynolds number gives access to 
smaller scales and shorter turn-over times with the result 
of having turbulence prevailing over waves and a tran- 
sition to a state that is consistent with a Kolmogorov- 
like scaling with dual energy and helicity cascade. Con- 
versely, at fixed Reynolds number, increasing rotation 
leads to a smaller inertial wave time scale and a preva- 
lence of rotation and of the e + h k, A regime found in 
Mininni and Pouquet (2009a) for helical rotating tur- 
bulence. A pure helicity cascade with no energy fiux 
is not observed, the highest value of the ratio of fluxes 
n_ff/[fc_FnB] being « 3.4. It is also remarkable that, when 
examining separately the energy and helicity spectra, in 
several cases we observe that the latter is well defined (ei- 
ther a Kolmogorov law or close to a k^^ law), while the 
index of the energy spectrum may not be so well defined. 

These results arc summarized in Fig. [6l which shows 
the runs that have e + /i ~ 4 with black dots, and the 
other runs with grey dots. The two axes correspond to 
the normalized ratio of the helicity to the energy fiux 
to small scales ni//[A:j;-n£;] as a function of the dimen- 
sionlcss parameter Nc — ReRo — U^^g/[iyH]. Since the 
effect of viscosity and rotation are antinomic, their prod- 
uct, properly adimensionalized by the kinetic energy of 
the system, is a determining parameter. Indeed, the 
clear demarcation discussed before can be observed as 
the abrupt partition of the black and grey dots around 
Nc ~ 500. This parameter is easily obtained in the 
case when the energy spectrum is E{k) (e£;r2)^/^fc~^ 
(DubruUe and Valdettaro 1992, Zhou 1995) by simply 
searching for the proportionality of the inverse of the dis- 
sipation wavenumber that emerges from the above spec- 
trum ko ~ [e£;/(i^^ri)]^/^ and the characteristic large 
scale of the flow (Canuto and Dubovikov 1997). This 
is consistent with the idea above that sufficiently small 
scales must be excited for the eddy turn-over time to be 
smaller than the inertial wave period and have a transi- 
tion towards the turbulent regime. Taking helicity into 
account, its flux en also comes into play and dimensional 
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FIG. 6: Scatter plot of the normalized ratio of helicity flux to 
energy flux as a function of Nc = ReRo = Urms/W^]- Note 
the rather sharp vertical partition around Nc ~ 500. 



analysis becomes undetermined, but the idea remains 
that the dissipation of energy evaluated using the appro- 
priate energy spectrum where waves play a role should 
lead to a parameter that governs the dynamics. 

This result may seem at odds with the conclusions 
drawn in Nolan (2005) where the importance of a dif- 
ferent combination of dimensionless numbers, namely 
Re/ Ro = flL^ /ly,,, is emphasized for studying the physics 
of tornadoes. This latter parameter can be viewed as 
the ratio of the vortex circulation to the eddy viscosity 
I/* and is sometimes called the vortex Reynolds number, 
balancing frictional and advective terms in the boundary 
layer. However, it should be noted that in our simulations 
there are no boundary layers (although internal bound- 
ary layers develop), and, furthermore, in the presence of 
strong turbulence, the eddy viscosity that replaces j^* in 
the above expression becomes proportional to U^^si ^^nd, 
thus, it may be that in the turbulent case the parame- 
ter ReRo = U^j^g/{vil.) becomes the relevant parame- 
ter. Moreover, note that the purpose of our parametric 
study is quite different: whereas Nolan (2005) seeks a 
predictive parameter before a tornado forms, we are here 
dealing with the overall scaling properties of small scale 
fluctuations in a rotating turbulent flow with updrafts or 
downdrafts (provided by the helicity). 

Besides different scaling laws, the other measurable dif- 
ferences we found between the two regimes identified in 
this paper, at fixed rotation rate, are a longer turn-over 
time tni^ in the case when e + /i = 4 together with a 
slightly slower growth rate in the inverse cascade, due 
presumably to a more efficient Beltramization of small 
scales because of the excess of small scale helicity. It 
would be of interest to investigate in detail the structure 
of strong laminar columns that form in these flows (Bel- 
trami Core Vortices, see Mininni and Pouquet 2009c), 
and relate them to their far-field environment, but this 



is left for future work. However, a preliminary examina- 
tion of these structures does not reveal any marked dif- 
ference of behavior between the two regimes: the same 
spatial juxtaposition of persistent laminar large colum- 
nar vortices elongated along the axis of rotation, and a 
web of intense small-scale vortices is visible in both cases 
(Mininni and Pouquet 2009bc). These structures are dif- 
ferent from what happens in a non-helical fiow (Mininni 
et al. 2008) in which case no laminar columnar vortices 
appear in the turbulent regime; this may indicate that, 
whatever the spectral law, helicity plays an important 
role in determining the statistical properties of rotating 
helical turbulence and the stability of its structures. 

Another remark that can be made at this point is that, 
as in all that precedes, the emphasis was on the impor- 
tance of helicity, which is a topological invariant measur- 
ing the knottedness of vortex lines (Moffatt 1969), but 
no computations were considered of rotating turbulence 
without net helicity. A continuation of this work will 
obviously involve performing computations with forcing 
possessing different degrees of relative helicity. In fact, 
in the non-helical case of a Taylor-Green forcing studied 
in Mininni et al. (2008), a spectrum was found, and 
one may ask whether it will also disappear in favor of a 
pure Kolmogorov spectrum as the Reynolds number is 
increased for fixed Rossby number. This point is left for 
future investigation. 



V. DISCUSSION AND CONCLUSION 

A significant effort has been put into modeling tur- 
bulent flows, both in the engineering and atmospheric 
contexts, as well as in astrophysics. What the results 
presented in this paper show is the fact that the he- 
lical model developed previously in Baerenzung et al. 
(2008a) works reasonably well when compared to high 
resolution high Reynolds number and moderate Rossby 
number direct numerical simulations of helical rotating 
turbulence. In this particular case, the inclusion of helic- 
ity improves the results, whereas in the non rotating case 
with both non-helical or helical forcing, or in the rotat- 
ing case with non-helical forcing, the inclusion of helicity 
in the model neither enhances nor degrades the results. 
For the same problem, the ChoUet-Lesieur (1981) model 
does not obtain the correct growth rate of energy in the 
inverse cascade; it can be viewed as somewhat dcflcient 
insofar as it seems too dissipative. This can be linked to 
the fact that, in the presence of strong helicity, the non- 
linearities are damped (Kraichnan and Panda 1988) and 
thus the turbulent dissipation is substantially diminished 
when compared to the non-helical case, particularly so in 
the rotating case (Teitelbaum and Mininni 2009). 

The CPU and memory usage savings when comput- 
ing with the model are impressive, since a grid of 96'^ 
points was used in all tests of the LES against the DNS 
run on a grid of 1536'^ points for the same Reynolds and 
Rossby numbers. It is possible that further gain may be 
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obtained when comparing to DNS performed on larger 
grids at higher Reynolds numbers, but such an evalua- 
tion of the optimum gain in a LES will have to wait for 
petascale computers and beyond. Indeed, a minimum of 
inertial range has to be resolved in the LES computa- 
tion in order to compute the eddy diffusivity and eddy 
noise expressed in terms of the resolved energy and he- 
licity spectra. Another advantage of the LES model is 
that doing temporal three-dimensional visualization of 
flows at resolution of 1024'^ grid points and above is still 
a costly exercise, demanding in memory usage as well as 
computer and human time. The LES runs, from that 
point of view, give a proxy visualization of the flow at 
almost no cost and are potentially very useful even when 
errors in some of the global statistical properties of the 
flow may be as large as 10%, as observed in some of the 
LES test runs on grids of 96^ points. 

A comparison with an under-resolved DNS on a grid 
of 160'^ points showed that all LES perform significantly 
better. Both the energy evolution and its spectral dis- 
tribution are wrong for the under-resolved DNS: there 
is insufhcient dissipation and the inverse energy cascade 
is too vigorous with an accumulation of energy at the 
largest scale as well as in the smallest scales, while the 
spectrum at intermediate scales is greatly under-valued. 
There is no doubt that at an equivalent resolution and/or 
numerical cost, the several LES tested here work signif- 
icantly better. Finally, we note that the LES-PH model 
has twice the cost in computational time compared to a 
DNS at the same grid resolution. 

An obvious application of the helical model tested here 
for rotating flows is to explore regimes unattainable to- 
day with DNS. We have done so in the second part of this 
paper, examining the scaling laws of energy and helicity 
in rotating flows in the direct cascade range. We find that 
there can be two regimes, one where inertial waves inter- 
act with helical eddies in the limit of low Rossby number, 
yielding energy spectra steeper than the Kolmogorov law, 
and one close to the Kolmogorov regime. 

Several questions need to be examined in the future, 
among which are: What are the structures in the regime 
that obtains at large Reynolds number? Are the prop- 
erties of fully developed turbulent flows at low Rossby 
number close to the predictions of weak turbulence the- 
ory (Galtier 2003)? How does the inverse cascade scale, 
and how does it saturate, when sufficient resolution is 
present at large scale? What will be the effect of adding a 
friction term at large scale? What is the effect of a mod- 
erate amount of helicity (only the cases of maximal or 
zero helicity forcing have been tested up to now)? Would 
a different choice of forcing (such as a two-dimensional 
force, or a random force) affect the results? 

Furthermore, the model presented in this paper does 
not include anisotropics or inhomogeneities of the small 
scales, nor does it take into account bottom topogra- 
phy (see e.g., Fredriksen 1999) or other realistic physics 
relevant when dealing with atmospheric flows. Memory 
effects are neglected as well in our modeling since the 



EDQNM on which it is based is a Markovianized clo- 
sure. In that light, a stochastic approach as developed 
by several authors (see e.g. Majda et al. 2003, Delsole 
2004, Fredriksen and Keppert 2006) has led to signifi- 
cant progress in the modeling of turbulent flows as tested 
against direct numerical simulations at moderate resolu- 
tions. Further improvements and tests will be needed 
to capture as well the memory effects of turbulent flows. 
Another obvious drawback of closure models of turbu- 
lence such as the EDQNM is that all information on mo- 
ments of the stochastic velocity field above second-order 
is lost, and phase information among Fourier modes is 
lost as well. Thus, for example, intermittency is not 
present in the EDQNM, although it is observed in the 
EDQNM-based LES since the LES, in principle, cap- 
tures sufficient information on the structure of the in- 
ertial range. Such improvements will require non-trivial 
developments. On the other hand, adding to the dy- 
namics one or several scalar fields, such as the potential 
temperature and the water vapor, cloud water and rain 
water mixing ratios as done in Wicker and Wilhelmson 
(1995), passively advected by the flow, is not necessar- 
ily that cumbersome, since the EDQNM for the passive 
scalar problem has been written and thus the transport 
coefficients are known. Using adaptive mesh refinement 
in the presence of boundaries, as done in Wicker and 
Wilhelmson (1995), or possibly with spectral accuracy 
(Rosenberg et al. 2006, 2007) will also enhance our ca- 
pacity to analyze complex flows. 

Some physical models have also incorporated helicity 
on the dynamics (Lautenschlager et al. 1988, Yokoi and 
Yoshizawa 1993, Li et al. 2006). For example, moti- 
vated by observations of tropical cyclones giving esti- 
mates for the averaged helicity in a variety of flows (An- 
thes 1982, Etling 1985, Lilly 1986), a hehcal subgrid- 
scale paramctrization was proposed in Lautenschlager et 
al. (1988), following similar studies in magnetohydrody- 
namics (see Krause and Riidiger 1974 for the case of neu- 
tral fluids). These analyses differ from the present study 
insofar as they concentrate on large-scale instabilities, 
whereas this paper has been devoted to the issue of small- 
scale statistics. For example, in Lautenschlager et al. 
(1988), the modeling, backed up by low-resolution DNS, 
writes in terms of transport coefficients proportional to 
the vorticity, inspired from the equivalent destabilizing 
effect of small-scale helicity on large-scale magnetic fields 
in MHD. Instabilities involving three derivatives of the 
velocity (in Fourier space, ~ fc'^{)(k, t)) in the presence 
of helicity were computed in Pouquet et al. (1978) us- 
ing the renormalization group when considering the limit 
fc — > 0, in which case they are sub-dominant. The prob- 
lem remains of properly modeling the effect of small-scale 
helicity on large-scale flows, the approach taken in this 
paper being a modification to the eddy viscosity integrat- 
ing a helical component, although large scale instabilities 
may develop as well (Frisch et al. 1984). 

Finally, our model is isotropic in the unresolved scales, 
an assumption that can of course be relaxed (see Cambon 
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and Scott 1999, Sagaut and Cambon 2008 and references 
therein). It would lead to more complex expressions for 
the energy transfer terms involving variations of corre- 
lation functions in terms of both k±_ and fcy (where _L 
and II refer to the direction of rotation). This implies 
that a numerical integration of such anisotropic closures 
is significantly more costly that what is performed in the 
work presented here, since angles as well as wavenumbers 
have to be discretized. As shown above, the isotropy as- 
sumption works relatively well at the moderate Rossby 
numbers tested here, which are close to atmospheric val- 
ues although at substantially lower values of the Rossby 



number, more complex models may have to be developed. 
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